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ABSTRACT 



We present a kinematic analysis of the globular cluster (GC) system in M31. 
Adding new velocity data for 150 GCs in our wide-field survey to those in the 
literature, we increase the number of M31 GCs with measured velocities by 42 per 
cent, to 504. Using the photometric and spectroscopic database of these 504 GCs, 
we have investigated the kinematics of the M31 GC system. We find that the all 
GC system shows strong rotation, with rotation amplitude of v rot ~ 190 km s _1 , 
and that a weak rotation persists even for the outermost samples at \Y\ > 5 kpc 
where Y represents the projected distance from the major axis. The rotation- 
corrected velocity dispersion for the GC system is estimated to be <7 pr ~ 130 km 
s _1 , and it increases from a p>r ~ 120 km s" 1 at \Y\ < 1 kpc to o V)T ~ 150 km 
s _1 at \Y\ > 5 kpc. These results are very similar to those for the metal-poor 



GCs. This shows that there is a dynamically hot halo in M31 that is rotating 
but primarily pressure-supported. We have identified 50 "friendless" GCs, and 
they appear to rotate around the major axis of M31, unlike M31's disk rotation 
around the minor axis. For the subsamples of metal-poor and metal-rich GCs, 
we have found that the metal-rich GCs are more centrally concentrated than the 
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metal-poor GCs, and both subsamples show strong rotation. For the subsamples 
of bright and faint GCs, it is found that the rotation for the faint GCs is stronger 
than that for the bright GCs. We have identified 56 GCs and GC candidates with 
X-ray detection including 39 genuine GCs with measured velocities. It is found 
that the majority of X-ray emitting GCs follow the disk rotation, and that the 
redder, more metal-rich, and brighter GCs are more likely to be detected as X-ray 
emitting GCs, as seen for GCs in early-type galaxies. We have derived a rotation 
curve of M31 using the GCs at \Y\ < 0.6 kpc, and it agrees well for the range 
of R = 20 — 45' to that based on other tracers except for the planetary nebulae. 
We have estimated the dynamical mass of M31 using 'Projected Mass Estimator 
(PME)' and 'Tracer Mass Estimator (TME)' as M PME = 5.5+q;3 x 1O U M out to 
a radius of ~ 55 kpc and M TME = 19.2^; 3 x 1O 11 M for a radius of ~ 100 kpc, 
respectively. We finally discuss the implication of these results and compare the 
kinematics of GCs with that of planetary nebulae in M31. 



Subject headings: galaxies: clusters: general -- galaxies: individual (M31) 
galaxies: kinematics and dynamics — galaxies: star clusters 



1. Introduction 



Globular clusters (GCs) are the oldest stellar systems in the universe and thus maintain 
a fossil record of the early history of galaxies. They are much brighter than stars and the 
spatial distribution of the GCs is much more extended than that of the halo stars in galaxies. 
Therefore they are an excellent tool to investigate the kinematics of the main components 
(stellar and dark matter halo, bulge, and disk) of galaxies and to estimate the mass of their 
host galaxies. 

M31 is a unique object as the nearest spiral galaxy in the Local Group, providing an 
excellent place to study GCs. It has the largest GC system of any galaxy wit hin severa l Mpc . 
There 
e.g. 



(e.g., ( 



ore, there are numerous studies of GCs in M31 since the first study bvlHubble I ( 19321) 



van den Bergh 1969; S argent et al. 1119 77: 



Battistini et al. 1987:lHuchra et al. 



Morrison et al. ll2004HBurstein et al. ll2004l:lHuxor et al. 



1991 



Barmbv et al. 



2005 



2000: 



Beasley et al. 



Perrett et al. 



2005 



2002 



Puzia et al. 



2005 



Gallcti ct al. 



2006; K im et al 1120071) and see 



extensive lists of earlier references in iHodge I ( 19921 ); Ivan den Bergh I ( 20001 )). In addition 



there were several surveys and studies of the M31 GCs with X-ray detection, fi nding that the 



X-ray emitting GCs are commonly associated with low-mass X-ray binaries ( Supper et al 



2001 



2004 



Kong et al. 



2002 



Pietsch et al. 



Di Stefano et al. 112002: Kaaret 



2005 



Fan et al. 



2005 



2002 



Williams et al. 



Trudolyubov et al. 



2005 



20061 ) 



20041 : iTrudolvubov fc Priedl 
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Kinematics is one of the notable features for the M31 GCs, and has long been a focus 
of study. It is known that the kinematics of the M31 GCs is much different from that of the 
Galactic GCs i n the sense that the former is dominated by rotation, while the latter shows 
weak rotation (IPerrett et al. Il2002l ) 



The first detailed kinematic analysis of t he M31 GC system was given by lHartwick fc Sargent 



( 11974J ) using radial velocities for 44 GCs in Ivan den Bergh I (119691 ) . They estimated the to- 
tal mass of M31, M tot = 3.4 ± 1.4 x 1O 11 M , and found that the metal- poor GCs have a 
larger velocity dispersion than the metal-rich GCs. iHuchra et al. I (119821 ) obtained digital, 
low-dispersion spectra of 61 GCs, and found that the GC system rotates with an amplitude 
of 160 ± 40 km s _1 . They also found a measurable but small metallicity gradient and no 
significant spectr oscopic difference betw een the GCs with X-ray detection and those without 
X-ray detection. iFederici et al. I (J199CH ) using the velocity data of 81 GCs obtained a total 
mass of 3.2 x lO n M for a ra dius of 16 kpc using a projected mass estimator. A comparison 
of the data with the model of iRichstone fc Tremaine I (11984J ) showed that the GCs have cir- 
cular o r isotropic orbits, and lead to a total mass of ~ 7 x 10 11 M Q based on circular orbits. 



Later, IFederici et al. I (119931 ) updated the total mass of M31 to be 5 — 8 x 10 M for a 



radius of 30 kpc using velocity data for 176 GCs. 

In the 1990's, metallicity information for a large sample of GCs (N > 100) became 
available, a nd it was possib l e to i nvestigate the kinematics of GC subsamples based on the 
metallicity. IHuchra et al. I (J19911 ) estimated the radial velocities and metallicities for 150 
GCs, and presented the detailed kinematics of the GC system using expanded subsamples 
of GCs based on metallicity. They identified a weak metallicity gradient as a function of 
projected radius. They found that, in the inner region (< 2 kpc), the metal-rich ([Fe/H]> 
—0.8) GCs exhibit a rapid rotation of 100 — 200 km s _1 , while the metal-poor GCs show 
no significant rotation. Howeve r, the rotation of the two subsamples could not be easily 
measured at larger radii. Later, IPerrett et al. I (120021 ) showed, using higher precision data 
of over 200 GCs ( typical velocity measurement error of 12 km s _1 ), that the entire M31 
GC system is highly rotation-supported. This is in strong contrast to the case of the Milky 
Way where the GC system composed mostly of halo GCs is pressure-supported, although 
a minor system of the bulge GCs shows some rotation. They also found that the centrally 
concentrated metal-rich GCs appear not to be flattened, but are analogous to the bulge 
population with a strong rotation , and t hat the less concentrated metal-poor GCs also show 
a rapid rotation. iGalleti et al. I (120061 ) obtained velocity data for 42 genuine GCs, and 
produced a large data set of 349 M31 GCs with up-to-date velocities by combining all 
information in the literature. 



Recently, large datasets of velocities for M31 GCs and several deep photometric studies 
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of the halo fields in M31 h ave provided important clues for understanding the evolution of 
M31. iBrown et al. I (120031 ) presented very deep photometry of stars in a halo field along 



the minor axis of M31 at a projected distance of 11 kpc (corresponding to a deprojected 
distance of 53 kpc) and found that the M31 halo includes a major intermediate-age metal- 
rich population and a significant old metal-poor population. They sugg ested from thi s result 



that M31 may have undergone a major merger 6 — 8 Gyr ago. Later IBrown et al. I (120071 ) 



suggested, from photometry of another minor-axis halo field at a projected distance of 21 
kpc, that the outer spheroid is dominated by the debris of earlier merging events that might 
have occurred more than 8 Gyr ago. 



Morrison et al. I (12004 ) suggested the exist ence of an old t hin d isk composed of GCs in 



M31, using the data for 321 GCs presented in iPerrett et al. I (120021 ) . Since the GCs in the 



thin disk have a metallicity distribution from [Fe/H] ~ —2.2 to above solar, and the metal- 
poor GCs are regarded as being as old as their Galactic counterparts, they concluded that 
the thin disk has remaine d und isturbed for ~ 10 Gyr or more. This is contradictory with the 

(200 3J) th at M31 may have undergone a major merger 6 — 8 Gyr ago. 



scenario of 



Brown et a 



Later iBeasley et al. I ((2004J, 120051 ) suggested that the majority of M31 thin disk GCs may 
be significantly younger than 10 Gyr, up to ~ 1 Gyr and metal-rich (—0.20 <[Fe/H]< 0.35). 
The existence of young and metal-rich GCs in the disk is also confirmed in other studies as 
well dBurstein et alTlbooi [puzia et al. Il2005h . 



On the other hand lAshman &: Bird I (119931 ) first tried to identify subgro ups in the M31 



GC sy stem using the positional and velocity information for 144 GCs in iHuchra et al 



(J1991I ). They suggested that the clustering of GCs in position and velocity space may be 
a surviving signatur e of gaseous clumps from which the galaxy halo of M31 formed. Later, 



Perrett et al. I (120031 ). using 301 GCs with measured velocity and metallicity in lPerrett et al. 



(120021 ). searched for the subclustering of M31 GCs in the parameter s pace of projected radius, 
radial velocity, and metallicity. Interestingly iMorrison et al. I (120031 ) found a possible tidally 
distorted satellite of M31, Andromeda VIII, from the identification of a group of 5 — 12 
planetar y nebulae (PNe), 1 — 3 GCs, and two HI cloud s. Recent kinematic studies of M31 



streams (jlbata et al. I l200ll ; iMcConnachie et al. 1 120041 ) started examin ing the GC 



including Andromed a VIII as a possible progenitor of the stream (e.g., iFont et al. Il2006al : 



groups 



Fardal et al. Il2006a[ ) 



Although there were numerous studies to obtain the radial velocities and the metallicities 
for M31 GCs, there are still many GC candidates for which spectroscopic information is 
needed. For example, there are 569 GC candidates for which r adial velocities are not available 
in a combined sample of M31 GCs in iGalleti et al. I (120061 ) . Recently we have performed 



a homogeneous photometric and spectroscopic, wide- field survey of GCs in M31 covering a 
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field of ~ 3° x 3° centered on M31. We presented the first result of our survey: a catalog of 
new GCs in M31 (JKim et al. 1 120071 . Paper I). This is the second in the series of papers on 
this survey. 

In this paper, we present the results of a kinematic study of the GC system in M31 
using the velocity data of 150 GCs newly obtained in our survey and 354 GCs previously 
known in the literature. Section 2 gives a brief description of the data used in this analysis, 
and section 3 investigates the kinematic properties for various subsamples of M31 GCs. We 
discuss primary results in §4. A summary is given in the final section. 



We adopted a distance to M31 of 780 kpc derived by iMcConnachie et al. I (120051 ) based 



on the method of the tip o f the red giant branch, which is similar to the value given by 
Freedman fc Madore I (119901 ). 770 kpc. At this distance one arcmin in the sky corresponds 
to 0.22 7 kpc. We adopted the center position of M31 as a = 00 h 42 m 44. s 3, S n = +41°16 / 09 " 
(J2000; ICrane et al. Ill992l ). Following the geometric transformation of lHuchra et al. I ( 19911 ). 
we define the X and Y coordinates, which are the projected distances along the major 
and minor axis, respectively. X and Y increase towards the northeast and the northwest, 
respectively, and are given by, 



x' = sin (a — a ) cos 5, 

sin S cos 5 — cos (a — a ) cos 5 sin So, 

X = x' sin (PA) + y' cos (PA), 

Y = -x' cos (PA) + y' sin (PA), 



(1) 
(2) 
(3) 
(4) 



where x' and y' are Cartesian coordinate s for a and S with a cen ter at («o, Sq), and PA= 37.7° 
is the position angle for the major axis (Ide Vaucouleurs I I1958T ) . 



2. Data 



2.1. Master catalog of M31 GCs 

First we used the data of GCs in M31 given in Paper I, which describes the details of 
our survey's photometric and spectroscopic observations, data reduction, and the data set. 
Here we only give a brief summary of the data set of M31 GCs. 

We selected GC candidates in Washington CMT\ CCD images covering a ~ 3° x 3° 
field centered on M31, obtained at the KPNO 0.9 m telescope. Spectroscopic observations 
were made using the Hydra multifiber bench spectrograph at the WIYN 3.5 m telescope 
in September 2000 and November 2001 for 748 GC candidates including 106 previously 
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known GCs. We determined the radial velocities for GC candidates by cross-correlating 
the candidate spectrum with those of two bright bonafide GCs in M31. Finally, we have 
independently found 1164 GCs and GC candidates, of which 559 are previously known GCs or 
GC candidates and 605 are newly found GC candidates. Among the GCs and GC candidates 
in our survey, there are 211 genuine GCs (class 1 in Paper I) with measured velocity. We 
have derived new velocity data for 150 GCs in our wide-field survey. 



Galleti et al.l (120041 ) compiled all known GCs and GC candidates in M31 in the literature, 



and pu blished the Revised Bologna Catalogue (RBC) with 1035 objects. Later. iGalleti et al. 



(120061 ) published an updated version of the Revised Bologna Catalogue V2.0 (RBC2) with 
1164 objects. In order to produce a master catalog of radial velocities for M31 genuine 
GCs, we have secured, from the RBC2, 349 genuine GCs (class 1 in RBC2) with measured 
velocities and 5 GC candidates (ID B366, G137, B102, B430, B134 in RBC2) with measured 
velocities that were classified as genuine GCs in our survey. There are 61 objects in common 
between this study and the RBC2, if we use a matching tolerance of 4". We have compared 
our velocity measurements with published velocity measurements (references are in RBC2) 
for these common objects, as displayed in Figure [U This shows that our measurements, 
with a typical error of 35.3 km/s, agree with the published values within the combined 
errors except for nine objects showing velocity differences of over 110 km s _1 . We present 
the identification in RBC2 for these GCs in Figure [H We tried to find any cause of the large 
difference for the nine objects, finding some explanation for five of them (B109, B124, B293, 
B315, and NB 21), but none for the other four(B347, B200, B331, and B387). The reasons 
for large velocity differences for five of these nine objects are as follows: 

1. B109 There are tw o velocity measurements for this objec t : —6 13 ± 24 km s _1 in 



Huchra et al. I (J199l|) and -372±12 km s" 1 in lPerrett et al. 1(120021) Since the velocity 



difference between two measurements is large and iPerrett et al. I (120021 ) gives higher 
precision data, RBC2 quoted —372 ± 12 km s -1 . However, the measu red velocity in 
our st udy is —568 ± 47 km s , which is consistent with the value in iHuchra et al 
Jl99lh . 



2. B124 There ar e two velocity measurements f or this object: —75 ± 22 km s 



m 



Jablonka et al. I ( 19981 ) and 70 ± 13 km s 1 in iBarmby et al. I ( 20001). Since t he ve- 



locity difference between these two measurements is large and IBarmby et al. I (120001 ) 
gives the value with the smaller error, RBC2 quoted 70 ± 13 km s -1 . However, the 
m easured velocity i n our study is —138 ±49 km s -1 , which is consistent with the value 



m 



Jablonka et al. I (I19981 ). 



3. B293 There are two velocity measurem ents for this ob j ect: —345 ± 35 km s l in 



Huchra et al. I (ll99ll ) and —467 km s l in lFederici et al. I (119931 ) . RBC2 presented the 
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value of a weighted mean (—424 ±23 km s x ), although the velocity difference is about 
f30 km s _1 . However, t he measured ve l ocity in our study is —534 ± 31 km s _1 , which 



is closer to the value of iFederici et al. I (119931 ) . 



4. B315 There are thr ee velocity measu rements for thi s obje ct: —291 ± 108 km s 1 in 



Huchra et al. I ( 



1991) 



|), -434 km s" 1 in lFederici et al. I Jl993J ). and -559 ±12 km s" 1 in 



Perrett et al. I (120021 ) . RBC2 pres ented the value of a weighted mean (—437 ± 18 km 



s l ) excluding the measurement in lPerrett et al. I (120021 ) . However, the measured veloc- 



i ty in our study is — 566±46 km s *, which is consistent with the value of lPerrett et al. 

(EooJ. 



5. NB21 The velocity m easurement in our study is —311 ± 35 km s 1 , while that in 



Jablonka et al. I (119981 ) is —773 ± 40 km s 1 . Since the pointing coordinate of our 



spectroscopic observation is different from the photometric center of this object by 
about 1.72 ", and this object is located in Field 25 in Paper I that is close to the M31 
center, the observed spectrum in this study might be significantly contaminated by the 
M31 stellar light. 

Since our measurements for B109, B124, B293, and B315 are similar to those of at 
least one value in the literature, we prefer to keep our measurements for those GCs at this 
time. We keep our measurement for NB21 for further analysis, since there is no convincing 
evidence to change our measurement at this moment. 

From the weighted linear fit discarding these nine objects with large velocity differences, 
we derive a transformation relation between the two systems, 



^This study = 1.09(±0.03) wrbc2 ± 37.8(±8.8) km s~\ rms = 37.8 km s _1 , N = 52. (5) 

The radial velocities for the GCs in RBC2 were transformed onto our velocity system 
using equation ([5]), and the transformed velocities were used for further analysis. We used 
our velocity measurements for the above nine objects. Finally, we produced a master catalog 
of radial velocities for 504 genuine GCs in M31. We derived 7\ magnitudes for 483 out of 
504 GCs in the m aster catalog from the CCD images used for our survey as described in 



Kim et al. I (120071 ). 



Figure 2 displays the errors of the measured velocities for the M31 GCs versus T\ 
magnitude (a) and versus R, a galactocentric radial distance corrected for the inclination of 
M31 using the equation R = [X 2 ± (Y/cos i) 2 ] 1 / 2 where i is the inclination angle, i = 77.7° 
(b). In Figure 2(a), the errors for 136 metal-poor GCs and 63 metal-rich GCs derived in 



this study are plotted with those for 174 metal-poor GCs and 58 metal-rich GCs derived in 
the previous studies. In Figure 2(b), the errors for 46 bright GCs and 99 faint GCs derived 
in this study are plotted with those for 155 bright GCs and 104 faint GCs derived in the 
previous studies (Division of the entire sample of M31 GCs into subclasses depending on 
the metallicity and magnitude will be described in detail in the following section). Mean 
errors for the measured velocities are 27 ± 18 km s _1 for all 504 GCs, 40 ± 16 km s _1 for 
211 GCs measured in this study, and 17 ± 12 km s _1 for 293 GCs measured in the previous 
studies. Mean errors for the measured velocities of the bright and faint GCs are 35 ± 12 
km s _1 and 40 ± 17 km s _1 for 46 bright GCs and 99 faint GCs, respectively, measured in 
this study, and 17 ± 13 km s _1 and 17 ± 9 km s _1 for 155 bright GCs and 104 faint GCs, 
respectively, measured in the previous studies. In the case of the previous measurements, a 
sig nificant fraction of th e measurements for faint GCs have typical errors of 12 km s" 1 given 



by iPerrett et al. I (120021 ). much smaller than the values i n other studies, which leads to a 



smaller mean value for faint GCs. Excluding the data by IPerrett et al. I (120021 ). the mean 



errors are 18 ± 15 km s _1 and 22 ± 13 km s _1 for bright GCs and faint GCs, respectively. 
Figure 2 shows that the mean errors for the measured velocities increase as the GCs get 
fainter, and that they are larger in the inner region where the background level is higher 
than that in the outer region of M31. 

The me a n val ues of the radial velocities are derived using the biweight location of 
Beers et al. I (Il990l ): Vp~ = —281 ± 13 km s _1 for 211 GCs measured in this study, and 



285 ± 8 km s for all 504 GCs in the master c atalog. These values are consistent with 



those in t he previous studies ( v p = — 284±9 km s 1 in lPerrett et al. Il2002l and v p = — 296±12 



km s * in lGalleti et al. 1120061 ) within the uncertainty, and are also compatible with the M31 
systemic velocity measured in this study (v p = — 290 ± 69 km s _1 ) and that in the literature 



—300 ±4 km s * in Ide Vaucouleurs et al. Ill99ll ) within the uncertainty. The larger 



error in our estimate of the M31 systemic velocity is primarily due to the fact that we used, 
as templates, the spectra of GCs with small velocity dispersion to estimate the velocity of 
the nucleus of M31 that has much larger velocity dispersion. 



Figure [3] shows the radial velocity histogram for all 504 GCs as well as 211 GCs measured 
in this study. It appears that the majority of all GCs are in the range —750 km s _1 < v p < 
+ 100 km s -1 , and the velocity distribution for all GCs is roughly symmetric with respect to 
the mean value of the radial velocities. Interestingly, there is one GC with v p < —850 km 
s _1 measured in this study (ID 106 in Table 5 of Paper I). This GC is in the northern disk of 
M31 (Field 13 in Fig. 1 of Paper I), and has Ti=18.94, C-Ti = 1.83, and v p = -890 ±38 km 
s _1 . Since the signal-to- noise ratio (S/N) of the spectrum for this GC is only ~5 (measured 
using H/3 and Fe5015 lines), a higher S/N spectrum is needed to verify its radial velocity 
and membership clearly. 
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The velocity dispersion derived using the biweight scale of iBeers et al. I (119901 ) is esti- 



mated to be a p = 185lg km s 1 for 211 GCs measured in this study, and a p 
for all 504 GCs in the master catalog. Both agree well within the errors. 



178+* km s" 1 



In Figure HI we show the spatial distribution of 504 GCs with measured velocities 
overlaid on a 4° x 4° optical image of M31 from the Digitized Sky Survey. The spatial 
distribution of 211 GCs measured in this study is globally similar to that of the 293 GCs 
measured in the literature. However, the GCs measured in this study are mainly located 
in the outer disk of M31 since the spectroscopic observation for the outer disk was done 
repeatedly (see Paper I for the configuration of the spectroscopic observations). 



2.2. Subsamples 

Metallicities of 199 GCs in M31 were derived from the line indices measured in the 
WIYN/Hydra spectra by Lee et al. (2007, in preparation, Paper III). We made subsamples 
based on the metallicities and T\ magnitudes for 504 GCs i n the master catalog f or further 



analysi s. We combined the metallicity data of 176 GCs in iBarmby et al. I (120001 ) and 194 



GCs in iPerrett et al. I (120021 ) with our data of 199 GCs (Paper III) in order to include the 



met allicity data in t he ma ster catalog. There are 32 objects in common bet ween this study 



andlBarmby et al. I (120001 ) . and 34 objects in common between this study and lPerrett et al 



( 20021 ). We have compared our metallicity measurements with those in these two previous 
studies for the common objects, as displayed in Figure It shows that our measurements 
agree well with the publis h ed on es overall, while the scatter in the comparison between this 
stu dy and IBarmby et al. I (120001 ) is smaller than that in the comparison between this study 
and IPerrett et al. I (120021 ) . We also derive t ransfo rmatio n relations between th e measurement 
in this study, and those in IBarmby et al. I (120001 ) and in IPerrett et al. I (120021 ) using the GCs 
with small measurement errors, 



[Fe/H] This study = 0.92(±0.02) [Fe/H] Barmby et al - 0.11(±0.02), rms = 0.19, N = 30, (6) 
[Fe/H] This study = 0.87(±0.03) [Fe/H] Perrett ct al - 0.12(±0.03), rms = 0.20, N = 25. (7) 



The measured metallicities for the GCs in IBarmby et al. I (120001 ) and IPerrett et al 



(120021 ) were transformed into our metallicity system using equations (jSJ) and (|7j), and the 
transformed metallicities were used for further analysis. We used our metallicity measure- 
ments for the common objects. Finally, we obtained metallicities for 431 GCs in total. The 
details of our metallicity data will be given in Paper III. 
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We present the met allicity distribut i on of M31 GCs in Figure Ufa). The KMM test 



based on the algorithm of I Ashman et al. I (119941 ) yields the result that the hypothesis that a 
unimodal distribution fits the data better rather than a bimodal distribution can be rejected 
at a confidence level of 99.9%. The bimodal test results in the following two groups: the 
metal-poor ([Fe/H]< -0.905) GCs (N = 310) with a mean value of [Fe/H]= -1.47 {a = 0.37) 
and the metal-rich ([Fe/H]> -0.905) GCs (N = 121) with a mean value of [Fe/H]= -0.62 
(a = 0.32). The KMM test with a trimodal distribution also rejects the hypothesis of a 
unimodal metallicity distribution at a confidence level of 99.8%. The trimodal test results 
in the following three groups: the metal-poor GCs (N = 183 GCs) with a mean value of 
[Fe/H]= —1.67 (a = 0.29), the intermediate metallicity GCs (N = 155 GCs) with a mean 
value of [Fe/H]= —1.14 (a = 0.26), and the metal-rich GCs (N = 93 GCs) with a mean value 
of [Fe/H]= —0.52 (er = 0.27). Although there is no reason to reject the trimodal metallicity 
distribution for M31 GCs (to be discussed in detail in Paper III), we use the two subsamples 
of metal-poor and metal-rich GCs from the bimodal distribution to increase the number of 
GCs per subsample for kinematic analysis in this study. 

We derived the mean values of the radial velocities: Up = — 294i xl km s _1 for 310 metal- 
poor GCs, and Up = — 2501*4 km s~ x for 121 metal-rich GCs. The value for the metal-poor 
GCs is consistent with that for the entire 504 GC sample (up~ = — 285l§ km s _1 ), but the 
value for the metal -rich GCs is significan tly smaller than that for the metal-poor GCs, which 



was seen earlier in lPerrett et al. I (120021 ) . 



For the velocity dispersion, we found that the dispersion for the metal-poor GCs (a p 



183^ 5 km s x ) is significantly larger tha n that of the metal- rich GCs (er. 



This is in contrast with the result of iPerrett et al. I (120021 ) that the dispersions for the 



1543 1 



km 



metal-poor and metal-rich GCs agree within the quoted errors (a p = 155 ± 7 km s x for 
the metal-poor GCs and a„ = 146 ± 12 km s _1 for the m etal-rich GCs) . We f ind that this 



difference may be due to the fact that the sample used by IPerrett et al. I (120021 ) was derived 
mainly from the disk of M31, while ours are based on a much wider area, including a larger 
fraction of probable halo GCs. 

The luminosity function for M31 GCs is shown in Figure [6] (b). We used only 483 GCs 
for which T\ magnitudes are available in our photometric catalog of M31 survey. Since the 
survey depth is not uniform across the galaxy due to varying background of stellar light, 
we use only the GCs brighter than T\ = 18.5 to make subsamples based on T\ magnitudes. 
Finally, we divide the GCs into 201 bright (T x < 16.9) and 203 faint (16.9 < T x < 18.5) GCs 
so that the number of GCs in each subsample is similar. 
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Results 



Using the master catalog of 504 GCs in M31, we have investigated the kinematic prop- 
erties of the M31 GC system for several cases of subsamples: all GCs, metal-poor and metal 
rich GCs, bright and faint GCs, and X-ray emitting GCs. Primary results of the kinematics 
of the M31 GC system derived in this study are summarized in Table 1. 



3.1. All GCs 

3.1.1. Kinematic Properties 

In Figured we present the spatial distribution of M31 GCs with the measured velocity. 
Here we adopted as the systemic velocity of M31 the value derived for all GCs in the previous 
section, —285 ± 8 km s _1 . Figure [7] shows clearly a rotation signature. The receding GCs, 
with velocities greater than the systemic velocity of M31, and the approaching GCs, with 
velocities less than the systemic velocity of M31, are spatially separated effectively by the 
minor axis, indicating a disk rotation. Interestingly, there are some GCs showing a retrograde 
rotation (filled circles at +X and open circles at —X) that will be discussed in §3.1.21 We 
present the radial velocities for M31 GCs versus the projected radii along the major axis (X) 
and along the minor axis (Y) in Figure [7| (b) and (c), respectively. A linear least-squares fit, 
passing through (X, v p — vmsi) = (0, 0), along the major axis results in v p — VM3i = 4.69+q;55X 
km s _1 and v p — VM3i = 0.31^'^y km s _1 , where X and Y are units of arcmin. The large 
value of the slope between v p — i>M3i and X indicates a strong rotation of the M31 GC system 
around the minor axis, while the small value of the slope between v p — Vm31 and Y indicates 
no significant rotation around the major axis. 

In Figure [BJ we present the radial velocities of GCs as a function of the projected 
distance along the major axis of M31 and the velocity histogram of the GCs for the samples 
of different distance bins along the minor axis in order to investigate any rotational variation 
along the minor axis. The mean radial velocity in a distance bin of 10' is overlaid by a large 
square with an errorbar representing the velocity dispersion. We fit the mean radial velocities 
with a straight line representing a solid-body rotation for the inner region. We select the 
appropriate boundary that gives the minimum \ value by moving the fitting boundary from 
10' to 50' symmetrically. Then, a flat rotation curve is connected to the solid-body rotation 
curve at that radius. The rotation amplitude that is inclination-corrected, is estimated by 
one-half of the velocity difference between the flat rotation curves. We also estimate the 
rotation amplitude using the velocity histogram for M31 GCs. We fit the observed velocity 
histogram using the sum of two or three Gaussian functions. Then, the rotation amplitude 
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is represented by one-half of the velocity difference between the two Gaussian peaks - one is 
the largest positive value and the other is the smallest negative value. We derive finally the 
rotation velocity corrected for inclination, v rot = v Pirot /sin i where v p>TO t represents a projected 
rotation velocity The rotation amplitude of all GCs in M31 estimated by fitting the rotation 
curve is v rot ~ 1881^8 km s _1 , while that from the velocity histogram is t> ro t ~ I88I33 km 
s _1 . These two values agree very well. The uncertainties on these (and following) estimates 
represent 68% (1 er) confidence intervals. We computed the uncertainties using 1000 artificial 
data sets constructed by randomly choosing GCs up to the number of GCs in the real data. 
The fitting procedure is performed on the 1000 trial data sets, the results are sorted, and 
the values corresponding to the 16th and 84th percentiles are identified. The uncertainties 
are defined as the offsets between these values and the values computed using the real data. 
We have also derived the rotation-corrected velocity dispersion of the all GC sample, that 
is, the dispersion about the rotation curve shown in Figure 8: a p>r = 134JJ km s _1 . 

The velocity histogram for the thin disk \Y\ < 1 kpc (d) shows not only two rotation 
peaks but also a velocity dispersion peak due to the central bulge of M31. The rotation 
amplitudes decrease as the distance along the minor axis of M31 increases. Interestingly, the 
rotation amplitude is not zero even for the outermost samples at \Y\ > 5 kpc, although the 
small size of the sample for outer GCs leads to a large uncertainty (Figure [8] (i) and (j)). 

This result indicates that the halo (or extended bulge) of M31 may be not totally 
pressure-supported but is rotating as well. This result is consistent with the idea that the 
M31 halo is dominated by a moderately rotating large bulge rather th an by a non-rotating 



halo far from the plane, which is based on the kinematics of PNe in M31 (jHurley-Keller et al. 



200 4J). However, th i s does not nece s sarily indicate that there is no classical halo in M31 (see 



Irwin et al. I (120051 ); iKalirai et al. I (120061 )). as discussed in the next section. 



On the other hand, the rotation-corrected velocity dispersion of the GCs increases from 
a P)T = 119lg km s _1 at \Y\ < 1 kpc to a p>r = 15llg km s _1 at \Y\ > 5 kpc. Thus the velocity 
dispersion becomes about twice as large as the nominal rotation amplitude at \Y\ > 5 kpc, 
showing that pressure support is much larger than rotation support in the outer halo. 



3.1.2. Friendless GCs 

Although the GC system in M31 is dominated by disk rotation, there are also seen some 
GCs that do not follow the disk rotation. To identify the GCs that do not follow the disk 
rotation, we sel ect 'friendless' GCs ad opting the method used for the study of the kinematics 



of M31 PNe in Merrett et al. I (120031 ) . 
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In summary, we select the friendless GCs that deviate from the mean velocity of neighbor 
GCs on the sky more than n x a, where a is the velocity dispersion of neighbor GCs. The 
number of GCs that defines the local environment is taken to be ~ [iV a n gcs(= 504)] 1 / 2 , and 
n = 2 is used in this study. Among all 504 GCs, we find 50 friendless GCs. The spatial 
distribution of the friendless GCs is shown in Figure [9j It is seen that they are located 
mostly in the disk. However, most of them are located in the quadrants of retro-grade 
rotation (the first and third) in the v p — VM3i versus X diagram. The friendless GCs in the 
second and fourth quadrants have mostly velocities larger than the largest velocities of the 
GCs following the major disk rotation. These friendless GCs may have different origin from 
the main population of GCs following disk rotation. 

To highlight the kinematics of the friendless GCs, we plot, in Figure [TD1 the spatial 
distribution of only the friendless GCs with the symbol size scaled according to their velocity 
deviations from the systemic velocity of M31. Interestingly, it appears that the friendless 
GCs rotate around the major axis of M31 unlike the rotation around the minor axis for the 
disk population. 



Merrett et al. I (120031 ) suggested that the friendless PNe in their sample may link to the 
Southern Stream and the Northern Spur since the positions and velocities of the friendless 
PNe are consistent with th e two streams. Recent l y, other stellar substructures in M31 were 



found by several authors (IFerguson et al. 1 120021 ; iMcConnachie et al. 1 12004 ; iFardal et al 



2006bl ). Since the kinematic information on stellar substructures in M31 is available, it 
would be interesting to compare the kinematics of the friendless GCs with those of stellar 
substructures, which may provide a clue to understand the motion and origin of the friendless 
GCs, but it is beyond the scope of this paper. 

In Figure [TTl we show the photometric properties of Xi magnitudes and (C — Ti) colors, 
and spectroscopic metallicities for the friendless GCs. The friendless GCs have a mean T% 
magnitude of < T\ >= 17.13 mag with a deviation of a^ = 1.47 mag, which is similar 
to the values of < T% >= 17.09 mag with a deviation of gt x = 1-30 mag for the normal 
GCs. The mean metallicities are also similar between the friendless and normal GCs, but 
the mean (C — Ti) color for the friendless GCs (< C — T\ >= 1.66 mag with (Jc-Tx = 0.42) 
is slightly redder than that for the normal GCs (< C — T\ >= 1.55 mag with oc~t x = 
0.60). We used a Kolmogorov-Smirnov (KS) test to determine whether the friendless and 
normal GC populations were drawn from the same distribution. The hypothesis that the two 
distributions are extracted from the same parent population can be rejected at the confidence 
level of ~ 86% for the T\ magnitude distribution, ~ 91% for (C — Ti) color distribution, and 
~ 1% for the metallicity distribution. Therefore, only Xi magnitudes and (C — Xi) colors 
appear to be slightly different. However, this problem needs to be investigated further using 
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the reddening-corrected magnitudes and colors. 



3.2. Metal-Poor and Metal-Rich GCs 



Figure [[2] shows the spatial and radial distributions of the 310 metal-poor and 121 
metal-rich GCs. The most striking difference in the radial distribution between the two 
subpopulations is that the metal-rich GCs show a much stronger central concentration than 
the metal-poor GCs, whi ch is consistent with the results based on the smaller data set in 
Huchra et al. I (119911 ) and lPerrett et al. I (120021 ). In particular the radial distribution of the 
metal-rich GCs show a very narrow peak with a half width of about 10 arcmin, showing that 
they belong to the central bulge of M31. Median values of the absolute distances along the 
major axis for the metal-poor and metal-rich GCs are 23.1' and 16.0', respectively, and those 
along the minor axis for the metal-poor and metal-rich GCs are 9.3' and 4.7', respectively. 
In addition, a small number of metal-rich GCs are found in the outer region, even at about 
70' from the major axis. There are 9 GCs within the boundary of the stan dard diameter 
in NGC 205. All of these GCs are metal-poor, as found by earlier studies (IPerrett et al 



20031 ). 



To highlight the difference in spatial distribution between the metal-poor and metal- 
rich GCs, we plot, in Figure [T31 the number density contours for the metallicity subsamples. 
Several features are noted in Figure [131 First, the spatial distribution of the metal-rich 
GCs shows a much stronger central concentration than that of the metal-poor GCs, as seen 
before. Second, the spatial distribution of the metal-rich GCs shows almost circular structure 
in the central region of M31, while that of the metal-poor GCs shows an extended structure 
elongated along the major axis with a similar ellipticity to that of the M31 disk. This 
difference is less clearly seen in the sample of bright GCs (T\ < 17). Third, there are several 
substructures in both subpopulations. The most notable metal-poor substructures are one 
in the region —40' < X < —20' and Y ~ —10', and one in the region —60' < X < —40' 
and Y « 0'. Thes e two correspond to the GC groups IDs 10 and 11, respectively, in 



Perrett et al. I (120031 ). These are close to, but not consistent with, the position of And VIII. 
One metal-rich GC and four metal-poor GCs are found within the standard diameter of And 
VI II. The substructure of metal-poor GCs in NGC 205 corresponds to the GC group ID 5 



in 



Perrett et al. I (120031 ) . Two notable substructures among the metal- rich GCs are on the 
major axis: one at X « —25' and the other at X « 20'. These correspond, respectiv ely, to 



the GC group ID 8 and the sum of IDs 2 and 3, respectively, in lPerrett et al. I (120031 ). One 



substructure containing both metal-poor GCs and metal-rich GCs at X ~ — 10' and Y gj 20 ' 



(seen only in the bright GCs) corresponds to the GC group ID 7 in IPerrett et al. I (120031 ). 
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Therefore most of the GC groups in iPerrett et al. I (120031 ) are confirmed in this study. 



fn Figure dH we present the radial velocities for the metal-poor and metal-rich GCs 
versus the projected radii along the major axis (X) and the minor axis (Y). A linear least- 
squares fit, passing through (X, v p —vm3i) = (0, 0), along the major axis results in v p —VM3i = 
+5.12±^X km s _1 for the metal-poor GCs and v p - v hn i = +5.28±\-jj£X km s _1 for the 
metal-rich GCs. A least-squares fit along the minor axis results in v p — VM3i — — O^O^'^Y 
km s _1 for the metal-poor GCs and v p — VM3i = +IA6-2.34Y km s" 1 for the metal-rich 
GCs. As known previously, both the metal-poor and metal-rich GCs show a strong and very 
similar rotation around the minor axis and no or a weak rotation around the major axis. 
The kinematic difference between the two subsamples is not significant at this stage. 

To investigate any kinematic difference between the two subsamples in detail, we plot, 
in Figure [13 the radial velocities for the metal-poor and metal-rich GCs as a function of the 
projected distance along the major axis of M31 and the velocity histograms. We also plot the 
mean radial velocities of the GCs in a distance bin of 10' along the major axis. We estimate 
the rotation amplitudes for the two subsamples using the same fitting procedure applied in 
Figure [SJ The rotation amplitude estimated by fitting the rotation curves (left panels) shows 
that the rotation of the metal-rich GCs (t> ro t = 191^37 km s" 1 ) is indistinguishab le from that 



of the metal-poor GCs (v Tot = 193+41 km s _1 ). However, IPerrett et al. I (120021 ) found that 
the rotation of the metal-rich GCs (i> rot = 160 ± 19 km s _1 ) is slightly stronger than that of 
the metal-poor GCs (t> rot = 131 ± 13 km s _1 ). 

It is worth noting that the rotation curve in the metal-poor GCs is asymmetric at 
\X\ > 30 arcmin. The metal-poor GCs at \X\ < 30 arcmin follow a solid-body rotation, 
while those at \X\ > 30 arcmin do not. Since the distribution of the metal-poor GCs is more 
extended than the metal-rich GCs, the asymmetric rotation curve is seen only for the metal- 
poor GCs. Contours in Figure [T5] show clearly this asymmetry at —60' < X < —30' where 
the mean radial velocities are larger than the opposite end (at 30' < X < 60'). Note also 
that there is a significant substructure at X « —25' with v p — i>M3i ~ —200 km s _1 , found in 
a wide-field survey. This corresponds to the most notable metal-poor substructure described 
above. Therefore these substructures are distinguishable not only in spatial distribution but 
also in kinematics. This substructure may be related with the Southern Strea m emanating 



from the southwest disk of M31 found in a wide field survey of red giants by llbata et al. 



(2001). 



In addition, we derived the rotation amplitude from the velocity histogram of the metal- 
poor GCs (panel (b)). The rotation amplitude of the metal-poor GCs is estimated to be 
v ro t = 15llg km s _1 , which is slightly smaller than, but agrees within the error to, the value 
derived from the v p — VM31 vs X diagram above. The velocity histogram of the metal-rich 
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GCs shows a single strong peak, showing that the sample is dominated by the GCs in the 
central bulge. So we did not try to derive a rotation velocity for the metal-rich GCs from the 
velocity histogram. We have derived the rotation-corrected velocity dispersions: o pr = 129lg 
km s _1 for the metal-poor GCs, and a PjT = 121JT 10 km s _1 for the metal-rich GCs. Thus the 
rotation-corrected velocity dispersion for the metal-poor GCs is very similar to that for the 
metal-rich GCs. 

In Figures [TBI and [T71 we also present the radial velocities of the metal-poor and metal- 
rich GCs as a function of the projected distance along the major axis of M31 and the velocity 
histogram of the GCs for the samples of different distance bins along the minor axis in order 
to investigate the rotational variation along the minor axis. The metal-poor GCs show that 
the rotation velocity is roughly decreasing with increasing distance from the major axis, but 
the uncertainties are too large to derive a reliable rotational velocity for the outer region at 
\Y\ > 3 kpc. The metal-rich GCs show a similar trend to that of the metal-poor GCs at 
\Y\ < 3 kpc, but they suffer more from the small number statistics problem. 

The rotation-corrected velocity dispersion of the metal-poor GCs increases from <j v<r = 
119^g km s" 1 at \Y\ < 1 kpc to cr p ^ = 145^} km s _1 at \Y\ > 5 kpc. Thus the velocity 
dispersion becomes about twice as large as the rotation amplitude at \Y\ > 5 kpc, showing 
that pressure support is much larger than rotation support in the outer halo of M31. 



3.3. Bright and Faint GCs 

In Figure [TBI we show the spatial distribution of the bright and faint GCs. It is seen 
that the bright GCs are more centrally concentrated than the faint GCs. Median values of 
the absolute distances along the major axis for the 201 bright and 203 faint GCs are 14.9' 
and 22.2', respectively, and those along the minor axis for the bright and faint GCs are 
7.1' and 8.0', respectively. Thus, the difference between the two subsamples in the spatial 
distribution along the major axis is significant, while that along the minor axis is not. This 
is mainly due to the fact that a significant fraction of bright GCs are located in the central 
bulge of M31 and that a large number of faint GCs in the central region of M31 were not 
detected in the current surveys. 

We present the radial velocities of the bright and faint GCs versus the projected radii 
along the major axis (X) and along the minor axis (Y) in Figure [T91 The linear least-squares 
fits, passing through (X, v p — VM3i) = (0,0), along the major axis result in v p — VM3i = 
+4,.98t° f 6 X km s" 1 for the bright GCs and v p - v M31 = +5.06±g;|gX km s' 1 for the faint 
GCs. The least-squares fit along the minor axis results in v p — f A/31 = -t-O.lSl^F km s _1 
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for the bright GCs and v p — VM31 = —0-79~^{{lY km s _1 for the faint GCs. Both the bright 
and faint GCs show a strong rotation around the minor axis and no or weak rotation around 
the major axis. The kinematic difference between the two subsamples is not significant. 

We also plot the radial velocities for the bright and faint GCs as a function of the 
projected distance along the major axis of M31 and the velocity histograms in Figure [20j 
We estimate the rotation amplitudes for the two subsamples using the same fitting procedure 
applied in Figure [HI The rotation amplitudes estimated by fitting the rotation curves (left 
panels) show that the rotation of the faint GCs (v Tot = 2091^ km s^ 1 ) is stronger than that 
of the bright GCs (v rot = I29I35 km s _1 ). In addition, the rotation amplitude from the 
velocity histogram shows again that the rotation of the faint GCs (t> rot 



2+23 



1781^ km s" 



is stronger than that of the bright GCs (v Tot = 144+ 4 g km s" 1 ). The difference in rotation 
amplitudes between the bright and faint GCs is closely related to the difference in their 
spatial distributions. Since the faint GCs are more extended along the major axis than the 
bright GCs, the rotation amplitude of the faint GCs is estimated to be larger than that of 
the bright GCs. It is noted that the rotation-corrected velocity dispersion for the bright 
GCs, a 



p,r 



146^7 km s 1 is much larger than that for the faint GCs, a p< 



p.r 



107+« km s" 



3.4. GCs with X-ray detection 



To investigate the kinematic properties of X-ray emitting GCs in M31, we cross-correlated 
the GCs and GC candidates from the optical survey with the X-ray sources in M31. We 
used as the op tical source list, the cata log of 1778 GC candi dates prepared combining our 



survey, RBC2. lHuxor et al. I (120051 ) and Mackev et al. I (120071 ). 



We made a master catalog of X-ray sources in M31 by combining; 560 sources in th e 
ROSAT survey (ISupper et al. Il200ll ). 204 sources in the Chand ra survey (jKong et al. 1120021 ). 
85 6 sources in the X MM-Newton survey (IPietsch et al. 1120051 ). and 62 X-ray emitting GCs 



in 



Fan et al. I (120051 ). If a GC candidate in the optical survey lies within a circle of 2<r 



positional uncertainty of an X-ray source, we regard it as a match. We used lcr positional 
uncertainty for the ROSAT sources due to its large positional uncertainty (5—48 arcsec of 
lcr with a median value of 7 arcsec) . In ad dition, we matched optical GC candida tes with 62 



X-ray emitting GCs in iFan et al. I (120051 ) using the identification of GCs, since iFan et al. 
(120051 ) presented only a catalog of X-ray emitting GCs. 



Finally, we identify 56 GCs and GC candidates (classes 1 and 2 in Paper I and RBC2) 
with X-ray detections (40 genuine GCs including 39 GCs with measured velocit ies, and 16 



GC candidates). There are 13 GCs that were identified as X-ray emitting GCs in lFan et al. 
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( 120051 ). but were not included in our combined optical GC catalog. Therefore, we increase 
the number of X-ray emitting GCs in M31 by 69, including 13 GCs that were not presented 
in our GC catalog. For further kinematic analysis, we use only the homogeneous data set of 
56 X-ray emitting GCs and GC candidates presented in our optical GC catalog. 



We show, in Figure ETJ the spatial distribution of 39 genuine GCs with X-ray detection 
and measured velocity in company with 1 GC and 16 GC candidates with X-ray detection but 
no measured velocity. It is difficult to make strong conclusions about the spatial distribution 
of the X-ray emitting GCs since the X-ray surveys did not uniformly cover the whole region 
in M31. However, it is clearly seen that a large number of the X-ray emitting GCs are 
located in the central bulge region of M31 and that they show a large velocity dispersion. 
In addition, it appears that the majority of X-ray emitting GCs follow the disk rotation of 
M31. 



In Figure 1221 we show the photometric properties of the X-ray emitting GCs compared 
with that of all 504 GCs. We used the KS test to determine whether the X-ray emitting 
and the entire GC population were drawn from the same distribution. Interestingly, the 
hypothesis that the two distributions are extracted from the same parent population can be 
rejected at the confidence level of ~ 99% for the distribution of Ti magnitude, [C — Xi) color, 
and the metallicity. This result based on the analysis of the expanded sample of M31 X-ray 



emitting GCs confirms the previous results (IBellazzini et al. I Il995l ; iDi Stefano et al. 1 12002 



Trudolyubov &: Priedhorsky 1120041 ; iFan et al. 1 120051 ) . The properties of X-ray emitting GCs 
in M31 are consistent with those of GCs in early- type galaxies: redder, more metal-rich, and 



brighter GCs are more likely to be de tected as X-ray emitting GCs ( iSarazin et al. 1 12003 



Kim et al. 112009 : ISivakoff et al. 1120061 ). 



3.5. Rotation Curve 



Rotation curves are useful for inferring th e amount and distribution of underlying dark 



matter in galaxies (e.g.. ISofue fc Rubin Il200ll ). For M31, there are several studies to d erive 



the rotation curve using various tracers: H II reg ions (Rubin fc Ford 



faraun I ll99ll: JCarignan et al. tod. I2OO7L CO jLoinard et al. Ill995h . PNe (JMerrett et al. 
2006J), and GCs JPerrett et al. Ibooj l. 



1970; Kent 



19891) . HI 



To obtain the rotation curve of M31 using the sample of 504 GCs in this study, we first 
selected 454 GCs following the M31 disk rotation by rejecting 50 friendless GCs defined in 
§3.1.21 Then, we selected 75 GCs with |Y/X| < 0.5 and \Y\ < 0.6 kpc that are dominated by 
the rotational motion, since the inclination of M31 to the line of sight is as high as i = 77.7°. 
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We display the rotation curve using the selected GCs in Figure 1231 We plot the rotation 
velocity corrected for inclination, v rot = {v p — VMs\){R/X)/sni i, where X is the projected 
distance along the major axis of M31 and R is the deprojected radius in the plane of the 
disk, R = [X 2 + (Y/cos 2) 2 ] 1 / 2 . For each GC, we use the nearest 21 GCs to estimate the mean 
velocity, taken as the rotation velocity at that radius, and to estimate the velocity dispersion. 
We stop the calculation when the number of the nearest GCs is less than 21 (R < 8.9' and 
R > 46.4'). Typical la uncertainties using the bootstrap procedure for the rotation velocity 
and for the dispersion are 12 km s _1 and 8 km s _1 , resp ectively. For comparis on, we derive 
also the rotation curve of PNe using the velocity data in iMerrett et al. I (120061 ) by applying 
the same procedure as above. For PNe, typical la uncertainties for the rotation velocity and 
for the dispersion are both 4 km s _1 . 

It appears that the rotation curve of the GCs in this study agrees well with those of 
other tracers in the rang e R = 20 / — 45'. However, the rotation curve of PNe derived 



using the velocity data of IMerrett et al. I (120061 ) is different from the other profiles in the 
similar radial range. The rotation velocity for PNe at R = 20 — 60' is on av erage about 50 



km s x smaller compared with those of other tracers, as seen in Figure 33 of IMerrett et al 



(120061 ) . They showed that the presence of significant random velocities at all radii account 
for the asymmetric drift between the rotational motion and the local circular speed for PNe. 
In the inner region (R < 20') of M31, the rotation curves are significantly different from 
each other. A careful study is needed for the rotation curves in the inner region due to the 
contribution of the non-circular motion of the bulge and the projection effect of an inclined 
M31 disk. 

The velocity dispersion profiles for the GCs and the PNe are similar in general except 
for the inner region at R < 15'. The velocity dispersion decreases from ~ 120 km s -1 at 
R = 10' to ~ 80 km s _1 at R = 30', and changes little in the outermost area. The velocity 
dispersion of the GCs seems to increase again at R ~ 45', while that of the PNe keeps slowly 
decreasing with R. More data on the GCs further out in the halo are needed to confirm this 
difference, which can probe whether a pressure-supported hot halo exists or not in M31. 



3.6. Mass of M31 



We have estimated the mass of M31 from the kinematic information of the GCs in this 
study using two methods: 'Projec ted Mass Estimator (PME) ' and 'Tracer Mass Estimator 
(TME)'. PME was introduced by lBahcall &: Tremaine I (119811) f or the case of test particles 



orbiting a point mass, and later was extended by iHeisler et all ( 119851 ) for those tracing the 
total mass. The mass based on PME is given by 
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M, 



PME 



NG 



N 



Yl v ^ R p^ 



(8) 



where N is the number of objects, G is the gravitational constant, R p is the projected 
galactocentric radius, and v is the line-of-sight velocity relative to the systemic velocity of 
the galaxy. The constant f p = 32/tt is used for the objects with an isotropic orbit, while 



f p = QA/n for the objects with a radial orbit. We obtain Mpme = 5.5_ a3 x 10 M Q using all 
504 GCs out to a radius of 242' ~ 55 kpc with an assumption of isotropic orbits. 



Later lEvans et al. I (120031 ) introduced a new method of TME for the case that the 
tracer population such as GCs and PNe does not follow the underl ying mass distr i butio n 



dominated by dark matter. They applied the TME to 321 GCs of iPerrett et al. I (120021 ). 
and obtained, for R p ~ 100 kpc, Mtme = 1-2 x 10 12 M Q that is composed of a rotation 
contributi on, M rnt = 3 x 10 n M^ and of a pressure contribution, M press = 9 x 1O 11 M . 
Similarly, iGalleti et al. I (120061 ) obtained a larger mass of Mtme = 2.4 x 1O 12 M using the 



TME with the expanded sample of 349 GCs. In this study, we first select 301 metal-poor GCs 
considered to be halo GCs, and choose a radial range where the surface number density of the 
halo GCs falls off as a power law. The surface number density of our metal-poor G Cs follows 
approximately a power law R~ 3 beyond 30', similar to that of lEvans et al. I (120031 ). Finally, 
146 halo GCs in the range of projected radius, R p 
We set a = (isothermal-li ke galaxy), 7 = 4, R 



30' - 152' (6.7 - 34.2 kpc) are selected. 
,,. m = 6.7 kpc, and -R PiOU t — 100 kpc using 
the same assumptions as in lEvans et al. I (120031 ). Then, we obtain for R p ~ 100 kpc, the 
total mass of M31, Mtme = 19.2^ 13 x 10 11 M Q that is composed of a rotation contribution, 
M rot = 6.2+^ x 1O U M , and of a pressure contribution, M press = 13.liJ;| x 1O U M . 



To compare our estimatio n with the results in 



cal mass estimates using GCs (IHartwick fc Sargent 



he literature, we plot the various dynami 



1974k Ivan den Bergh Ill98ll; iFederici et al 



1993; Evans fc Wi 



kinson 



20001: 



Perrett et al. 



satell ite galaxies dCourteau fc van den Bergh 



20001 ). and streams (IFardal et al. Il2006ai) of M31 in Figure 



of the mass model for M31 are plotted (IWidrow et al. 112003 ; 



2002; 



199 



9: 



vvans et al. 



2003; 



Galleti et al. 



Evans fc Wilkinson 



200 



0; 



2006h. 
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241 In addition, several profiles 



Ibata et al. 112004 ; ISeigar et al 



20061 ; iGeehan et al. 1 120061 ) . The mass est imate based on PME derived in this study agrees 



well with other observati onal estimates by lEvans fc Wilkinson I (120001 ); IGalleti et al. 
and two mass models by lWidrow et al. I (120031 ); ISeigar et al. I (120061 ). 



(l2006h 



However, the mass estimates based on TME derived in this study and in IGalleti et al. 



( 120061 ) are larger than other estimates, and they are also larger than any models. If the TME 
is correct , the mass model of M3 1 needs a revision to a larger mass. This is also consistent 

(l20071 )'s suggestion based on the study of a new dwarf galaxy, 



with the 



Majewski et al. 
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AndXIV. Noting that AndXIV shows a large velocity of —203 km s 1 relative to M31 at a 



projected galactocentric distance of 162 kpc, iMajewski et al. I (120071 ) suggested that recent 



mass models of M31 need a revision to larger, if AndXIV is bound to M31. Similar discussion 
on the increase of the Milky Way Galaxy's mass base d on the mo s t dist ant satellite (d = 270 
kpc), the Leo I dwarf spheroidal galaxy, was given in iLee et al. I (119931 ). 



4. Discussion 



4.1. A dominant bulge versus a wimpy halo in M31? 



The existence of a halo in M31 has been controversial. While iHuchra et al. I (Il99ll ) 



suggested that the metal-poor GCs in M31 show little evidence of rotation, indicating that 
they belong to the halo of M31, other studies based on surface photometry, counts of red 
giant branch stars, and counts of PNe along the minor axis of M31 found an extended 



bulge, but no evidence of any halo component in M31 (IPritchet fc van den Bergh 



Mould & Kristian Ill986l : iHurlev-Keller et al. 112004 ; llrwin et al. 



2005 



Merrett et al. 



1994 



2006|). 



Howeve r, recent studies covering a very wide field of M31 are finding evidence for a huge halo 



in M31 (llrwin et al. 



2006; Gilbert et al. 



2005 



40 kpc from the center of M31 (llbata et al. 1120051 ) 



Guhathakurta et al. 112005k IChapman et al. 1120061 ; iKalirai et al. 



20061 ; llbata et al. 1120071) as w ell as a vast extended disk reaching about 



Irwin et al. 



(120051 ) presented a minor-axis surface brightness profile reaching out to a 
projected radius of 4° (« 55 kpc), combining surface photometry and faint red giant star 
counts. They found that the minor-axis profile can be fit by a de Vaucouleurs law with 
effective bulge radius R e s = 0°.10 (1.4 kpc) out to a projected radius of 1°.4 (« 20 kpc), 
and that there is a separate component beyond 1°.4 that can be fit by either a power law of 
index about -2.3 or an exponential law of scale length 14 kpc. One thing interesting about 
this component is that it has the same color as the inner region, indicating that it has a 
similar metallicity to that of the bulge (assuming the two components have the same age). 
In addition, the photometry of remote fields along the minor axis of M31 showed that most 
stars in these fields have unexpectedly high metallicity, [m/H]~ —0.5 ([Fe/H]~ —0.8) while 



there are some metal-poor stars as well ( Mould fc Kristian II 19861 ; iDurrell. Harris, fc Pritchet 



20011 ). Thus there are two components in the outer region of M31 that are distinguished in 
structure, but not in metallicity. We may say from the results based on stars that M31 has 
a huge dominant bulge reaching out to 4°, or that M31 has two bulges: an inner bulge (out 
to 1°.4) and an outer bulge (out to 4°). Even the inner bulge is much larger than the central 
bulge in our Galaxy, and the outer bulge may be called a halo. 
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The progenitors of PNe are low to intermediate mass stars so that they are often con- 
sidered as a good tracer for old stellar populations. For example, the radial number density 
profiles of PNe in M31 show excellent agreement with the - R-band surface brigh tness profiles 
of the stellar bulge and disk (IWalterbos fc Kennicutt II 198 71 ; iMerrett et al. 1120061 ) . In particu- 
lar, the number density pr ofile of PNe along the minor axis of M31 is consist e nt with the R 1 ^ 



profile out to one degree (IMerrett et al. 1 120061 ; iHurley-Keller et al. 1 12004 ) . IMerrett et al. 



( 120061 ) pointed out that there is little excess of PNe over the R l l 4 profile even at ~10 R e ff 
(one degree), indicating no separate halo component. This is also consistent with the re- 



sult based on the surface brightness profile along the minor axis (IPritchet fc van den Bergh 



1994 ; llrwin et al. 1120051 ). Considering the rotating kinemat i cs of the outer PNe at \Y\ > 6 
kpc in M31 as well as the above points, IHurley-Keller et al. I (120041 ) suggested that the bulge 
is the most dominating component in M31, while a nonrotating, pressure-supported halo is, 
if any, only a very minor component of M31. 

In general, or at least in our Galaxy, metal-poor GCs are considered as a tracer of the 
halo, while metal-rich GCs belong to t he bulge or t hick disk like the PNe, although the 
bulge does have some metal-poor GCs (JHarris II1996J ). Note that the kinematics of 17 old 
star clusters in M33, irrespective o f their metallicity, show evidence for the existence of a hot 
halo in M33 (jChandar et al. II2002I ). Therefore we may use metal-poor GCs to investigate the 
properties of any halo in M31. The kinematics of the metal-poor GCs along the major axis at 
\Y | < 1 kpc show a rotation as strong as the metal-rich GCs. The rotation amplitude of the 
metal-poor GCs decreases with increasing distance from the major axis, and the metal-poor 
GCs show some rotation (60 ~ 100 km s" 1 ) even far from th e major axis a t \Y\ > 5 kpc. 



These results are consistent with those of iPerrett et al. I (120021 ). but not with lHuchra et al. 
jl99lh . 



However, the rotation-corrected velocity dispersion for the metal-poor GCs is about 
twice as large as the rotation amplitude in the outer region at |Y"| > 5 kpc. The rotation- 
corrected velocity dispersion for the metal-poor GCs at \Y\ > 5 kpc in M31, a v ^ r = 1451}^ 
km s _1 , is comparable to that of the meta l-poor GCs in the Galactic halo at Rgc > 4 
kpc, o = 123 km s _1 given by ICote I (119991 ) . Considering that the spatial distribution of 
the metal-poor GCs in M31 is much more extended than that of the metal-rich GCs and 
that the velocity dispersion of the metal-poor GCs in the outer region of M31 is about 
two times larger than the rotation amplitude, we conclude that there is a dynamically hot 
halo in M31 that is rotating but primarily pressure-supported. This is consistent with the 
discovery of a pressure-supported metal-poor stellar halo based on the observations of red 



giant branch stars at rad ii from 10 kpc to 165 kp c (jChapman et al. 1 120061 ; I Gilbert et al. 



2006 



Kalirai et al. 1120061 ). I Chapman et al. I ( 20061 ) found from the analysis of a large sample 
of red giant branch stars: (a) there is seen a stellar halo component with low metallicity ( 
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[Fe/H]~ —1.4 ) at radii from 10 to 70 kpc, (b) there is little radial gradient of metallicity, 
(c) there is no significant evidence for rotation, and (d) the velocity dispersion of this metal- 
poor halo can be represented by cr v (R p ) = 152 — 0.90i? p km s _1 where R p is a projected 
galactocentric distance in units of kpc. 

In Figure 25 we display the radial variation of the mean metallic ity for the metal-poo r 
GCs in comparison with that for the red giant branch stars given by iKalirai et al. I (120061 ). 
A few features are noted in Figure 25: (a) the mean metallicity for the metal-poor GCs 
is much lower than that for the red giant branch stars at the same distance from the M31 
center; (b) the mean metallicity for the metal-poor GCs in the inner region is similar to 
that for the red giant branch star s much further out a t R p ~ 80 kpc that was considered to 
be the part of the stellar halo by IKalirai et al. I ( 20061 ); and (c) the mean metallicity of the 
metal-poor GCs shows little radial dependence: [Fe/H]= —1.48 ± 0.27 at R = 7.1 kpc to 
[Fe/H]= -1.62 ± 0.25 at R = 25.5 kpc (if we use all GCs, we obtain [Fe/H]= -1.24 ± 0.27 
at R = 6.7 kpc to [Fe/H]= —1.62 ± 0.25 at R = 25.5 kpc). These show that the metal-poor 
GCs are tracing the metal-poor halo, even in the inner region of M31. 



Recently llbata et al. I (120071 ) found a faint, smooth and extremely extended halo compo 



nent , reaching out to 150 kpc, about three times larger than the limit studied by llrwin et al 
( 120051 ). In addition, they found, from the photometry of red giant branch st ars, that the 



major ity of the stars in this extended halo are metal-poor, in contrast to the llrwin et al 



( 120051 ) finding that its integrated color is similar to the inner bulge. In spite of its huge size, 
the total luminosity of this halo derived assuming that it is symmetric is found to be ~ 10 9 
L Q , similar to the stellar halo of our Galaxy. Therefore it will be interesting to find GCs and 
investigate their kinematics in this outer halo of M31. 



4.2. Comparison of GC and PN kinematics 

PNe are an excellent tool for studying the kinematics of the bul ge and halo of a galaxy 



Recently huge data s ets for the PNe in M31 were published by IMerrett et al. I (120061 ): 



Halliday et al. I (120061 ). How do the kinematics of the GCs compare with that of the PNe 
in M31? We h ave derived the kinem atics of about 2500 PNe using the same method to the 
data given by IMerrett et al. I (120061 ). as used for the analysis of the GC data, which are 
summarized in Table 2. 

Figure 26 displays radial velocities (dots) versus the projected distances along the ma- 
jor axis (X) a n d the velocity histogram for the PNe in M31 based on the data given by 



Merrett et al. I (120061 ). like Figure 8 for GCs. We plot also the mean radial velocity and its 
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dispersion of the PNe in a distance bin of 10 ' along the major axis. Velocity histograms of 
PNe are displayed with Gaussian fits as in Figure 8. Figure 26 shows that the rotational 
properties of the PNe are very similar to those of the GCs. The rotational amplitudes of all 
PNe are v TOt ~ 1871*3 km s _1 from the rotation curve, and v rot ~ 198lg3 km s _1 from the 
velocity histogram, which are similar to those of GCs derived in this study, v rot ~ 188^28 km 
s _1 , and v TOt ~ 188i|§ km s _1 , respectively. The rotation amplitudes of the PNe decrease as 
the distance along the minor axis of M31 increases, but the rotation amplitude is not zero 
even for the outermost samples at \Y\ > 5 kpc, which are consistent with the results based 
on GCs derived in this study. 

The rotation-corrected velocity dispersion for all PNe is estimated to be a PtT = 94lJ 
km s _1 , which is much smaller than that for all GCs, a pr ~ 134 ± 5 km s _1 . The rotation- 
corrected velocity dispersion for the PNe at \Y\ < 1 kpc is o p ^ r = 971^ km s _1 , and it 
becomes smaller and shows little variation with \Y\: a p>r = 751^ km s _1 for 1 < \Y\ < 3 
kpc, a Ptr = 88l 4 km s" 1 for 3 < \Y\ < 5 kpc, and a Pjr = 79^ km s" 1 for \Y\ > 5 kpc. 
Thus the rotation-corrected velocity dispersion for the PNe at \Y\ < 1 kpc is consistent with 
that for the metal-rich GCs at \Y\ < 1 kpc, <y pr = 107^22 km s _1 . We cannot compare the 
kinematics of the PNe at \Y\ > 1 kpc with that of the metal- rich GCs at \Y\ > 1, because 
there are too few such GCs in this region. 

The rotation-corrected velocity dispersions for the PNe at \Y\ > 1 kpc are significantly 
different from those for the metal-poor GCs. That is, the rotation-corrected velocity disper- 
sion for the PNe at \Y\ > 1 kpc are smaller than that for the metal-poor GCs at \Y\ > 1 
kpc, and the former shows little variation depending on \Y\ for \Y\ > 1 kpc, while the latter 
increases with \Y\. These differences show that the PNe are tracing well the bulge, while 
the metal-poor GCs are tracing better the metal-poor halo, although their spatial coverages 
in the sky are similar. Therefore we conclude again that there do exist both an extended 
rotating bulge and a dynamically hot metal-poor halo in M31. 



4.3. Formation and Evolution of M31 

Main features related with the kinematics of the M31 GC system revealed in this study 
can be summarized as follows. Most of the GCs show a strong rotation around the minor 
axis. The rotation velocity decreases with increasing distance from the major axis. It is 
still as large as 60 ~ 100 km s~ x even at \Y\ > 5 kpc, indicating the bulge of M31 is 
large and rotation-d o mina ted, which is consist e nt wi th the finding based on the PNe by 



Hurley-Keller et al. I (12004 ) and iMerrett et al. I (120061 ) . The rotation velocity of the metal- 



rich GCs is similar to that of the metal-poor GCs in the region close to the major disk. The 
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spatial distribution of metal-rich GCs shows a stronger central concentration than that of 
the metal-poor GCs. The rotation-corrected velocity dispersion of the metal-poor GCs in 
the outer region of M31 is about twice larger than the rotation amplitude. There are several 
substructures in the spatial distribution of both metal-poor GCs and metal-rich GCs. There 
are about 50 friendless GCs all over the disk of M31, and most of them show retrograde 
motion. Interestingly it appears that these GCs show a rotation around the major axis, 
while most of M31 GCs show a rotation around the minor axis. 

What do these results in addition to the known properties of M31 tell us about the 
formation and evolution of M31 in comparison with our Galaxy? The most distinguishable 
kinematic difference between the M31 GCs and the Galactic GCs is a difference in the bulge 
and halo kinematics. In our Galaxy metal-poor halo GCs are the most dominant population 
among the GCs, and metal-rich GCs in the bulge (or thick disk) are a minor population . 
The bulge GCs show significant rotation (v ro t = 193 ± 29 km s" 1 ) , while the halo GCs are 
press ure-supported (v rat = 43 ± 29 km s _1 and a v = 123 km s _1 ) (jArmandroff Ill989l ; ICote 



19991 ). In contrast, the disk and bulge GCs are the most dominant population in M31, and 
the halo GCs may be a minor population compared with the Galactic halo GCs. 

The current picture of the formation and evolution of the halo in our Galaxy involves 
basically two classical scenarios: (a) The monolithic collapse model according to which 
the stellar halo is formed by the rapi d collapse of a monolithic protogalact ic cloud in a 
dynamical time scale of about 0.1 Gyr (jEggen. Lynden-Bell. fc Sandage Ill962l ). and (b) the 
accretion/merging model according to whi ch the stellar ha l o is formed via successive ac cretion 
and merging of protogalactic fragments (ISearle fc Zinn I Il978l ; IWhite fc Rees 1 119781 ). The 
latter is related with the current paradigm for explaining the formation of the cosmological 
large scale structure: a hierarchical merging scenario in the ACD M universe. It appear s 
that both scenarios are involved in forming the halo of our Galaxy (IChiba fc Beers 1120001 ). 
While the inner part of our Galaxy can be explained mainly by the first, the outer part of 
our Galaxy is explained better by the latter. 



Recently, iRenda et al. I (J2005af ) suggested from a semi-analytic multizone chemical evo- 
lution model that the observed higher metallicity in the M31 halo compared with the Galactic 
halo can be explained by two scenarios: (a) a higher star formation efficiency for the halo 
that is expected from the accretion of smaller satellites, or (b) a larger reservoir of infalling 
halo gas with a longer halo phase, producing the intermediate-age populatio n in the halo. 



These w ere also confirmed in the semi-cosmological spiral galaxy simulation by lRenda et al. 



( J2005bl ). indicating that the M31 halo may have a more protracted as s embly history com 



pared with the halo of our Galaxy. On the other hand, iFont et al. I ( J2006bl ) found from 
the simulation of the stellar halo in our Galaxy based on the hierarchical scenario: (a) the 
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stellar halo of our Galaxy formed inside out via accretion of satellites; (b) the inner halo at 
R < 20 kpc assembled rapidly, and most of its mass were in place more than 8-9 Gyr ago; (c) 
satellites surviving today were accreted within the past few Gyr; and (d) a major fraction of 
the present-day stellar halo was originally from a few massive satellites with M ~ 10 8 — 10 10 
Mq. They also pointed out that the much lower metallicity of the Galactic halo compared 
with the M31 halo indicates that our Galaxy could have formed mainly through smaller 
satellites. Noting that our Galaxy has a stellar mass a nd angular momentum that are 2 to 3 
times lower than M31, which is a typical spiral galaxy, iHammer et al. I (120071 ) suggested that 
our Galaxy has an exceptionally quiet formation history involving no major merger during 
the last 10-11 Gyr. However, it is noted that the term 'M31 halo' used in these studies cor- 
responds to the extended bulge, not to the metal-poor halo in M31. Therefore these studies 
explain the formation of the extended bulge in M31 and the halo in our Galaxy. 

M31 might have formed and evolved in a way similar to our Galaxy, with some dif- 
ferences. The existence of the rotating metal-rich central bulge, the radial gradient of the 
metallicity and the extended smooth metal-poor halo in M31 can be explained by a rapid 
dissipative collapse of a proto- Galactic cloud, while the extended bulge and i ts several sub 



struct u res including giant streams and friendless GCs (jlbata et al. I (120051 ); I Gilbert et al. 



(120071 ) ; llbata et al. I (120071 ) and this study) , are evidence for accretion and merging of small 
scale fragments. Differences between M31 and our Galaxy are that the degree of merging 
was larger, the scale of merging fragments was larger, and the duration of merging (or the 
frequency of merging) was larger in M31 compared with our Galaxy. Remaining questions 
are: (a) Why does M31 have an extended metal-rich bulge in M31, which is not seen in our 
Galaxy?; (b) Why is the M31 extended bulge a strong rotator?; (c) Are there metal-poor 
GCs that do not show systematic rotation in the halo of M31?; (d) Why do the metal-poor 
GCs in M31 show some rotation while those in our G alaxy do not?; and (e) Is there indeed 
a population of young disk GCs in M31 as posited by iMorrison et al. I (120041 ). which has no 
counterpart in the Galaxy? 



5. Summary 



We produced a master catalog of 504 G Cs with measured velocity in M31 combining 
the data for 211 GCs measured in our survey (IKim et al. 1120071 ) and 354 GCs (including 61 
common GCs) available in the literature. Using the photometric and spectroscopic database 
of these GCs, we have investigated the kinematics of the GC system of M31. Our velocities 
are in good agreement with published values for objects in common, with a typical error of 
35.3 km/s. Primary results are summarized below: 
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1. The mean value and dispersion of the radial velocities for the 211 GCs measured in 
this study are Up~ = — 281^}g km s _1 and o p = 185^g km s _1 , respectively. Using all 
504 GCs in the master catalog, we obtain the mean value and dispersion of the radial 
velocities as Vp~ = — 285+g km s _1 and a p = 178JI 4 km s _1 , respectively 

2. For all GCs, the rotation amplitude estimated by fitting the rotation curve is v TOt = 
I88I28 k- m s~\ while the one from the velocity histogram is v rot = I88+33 ^- m s_1 - The 
rotation amplitude decreases as the distance along the minor axis of M31 increases. 
However, the rotation amplitude is not zero for the outermost samples at |V| > 5 kpc, 
indicating an extended rotating bulge or a hot rotating halo. 

3. We have identified 50 friendless GCs among all 504 GCs. It appears that the friendless 
GCs rotate around the major axis of M31 unlike the rotation around the minor axis 
for the disk population, which calls for more sophisticated investigation. A KS test 
yields that the difference between the friendless and normal GCs is significant for the 
distribution of T\ magnitude and (C — T\) color. 

4. For the subsamples of the metal-poor GCs and metal-rich GCs, it is found that the 
metal-rich GCs are more centrally concentrated than the metal-poor GCs, and that the 
rotation amplitudes for the metal-poor GCs and metal-rich GCs are similar. For the 
subsamples of bright and faint GCs, it is found that the bright GCs are more centrally 
concentrated than the faint GCs, and that the rotation for the faint GCs is stronger 
than that for bright GCs. 

5. The rotation-corrected velocity dispersion for all GCs is estimated to be a PtT ~ 130 
km s _1 , and it increases from o v , r ~ 120 km s _1 at \Y\ < 1 kpc to o p , r ~ 150 km 
s" 1 at \Y\ > 5 kpc. These results are very similar to those for the metal-poor GCs. 
This shows that there is a dynamically hot halo in M31 that is rotating but primarily 
pressure-supported. 

6. We have identified 56 GCs and GC candidates with X-ray detection including 39 gen- 
uine GCs with measured velocities. Kinematic difference between the X-ray emitting 
and other GCs is not significant at this stage, but further studies are needed with an 
expanded sample of X-ray emitting GCs. The photometric properties of the X-ray 
emitting GCs show that the GCs that are redder, more metal-rich, and brighter are 
more likely to be detected as X-ray emitting GCs, as seen for GCs in early-type galaxies 

7. We derived the rotation curve of M31 using the metal-poor GCs at \Y\ < 0.6 kpc 
from 8.9 to 46.4 arcmin. The rotation curve using the GCs in this study is consistent 
with those using othe r tracers in the range 20 — 45' except for that based on PNe of 



Merrett et al. I fj2006h . 
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The dynamical mass of M31 using the kinematic information of the GCs in this study 
is estimated to be M PME = 5.51q; 3 x 1O 11 M out to a radius of R ~ 55 kpc using the 
PME, and M TME = ^-Z-H x 1O U M at R ~ 100 kpc using the TME. 
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Table 1. Kinematics of the M31 Globular Cluster System 



1*1 N V a b p a c p>r < M v^ 2 

(kpc) (km s _1 ) (km s _1 ) (km s _1 ) (km s _1 ) (km s -1 ) 

All GCs 

~~0< 504 -285±g 178±| 134+^ 188+g 188l| 

0-1 149 -287lljj 200+^ 119+° 188+^ 234+^ 

1-3 203 -29711? 172l| HOljj 180l|j 1471}? 

3-5 65 -2861H 156±i| 1411U 1221^ 89l^ 

5< 87 -257111 1621H 151±| 98l^? 59lg 

Metal-Poor GCs ([Fe/H]< -0.905) 

~~0< 310 -2941H 1833 129+J 193lJJ 151±| 

0-1 79 -2961^ 210l| 119ll£ 297llj£ 239±J| 

1-3 135 -306111 175l£ lllllg 2061^ 152±1| 

3-5 37 -30213° 1671^ 20ll^ 30111$ 64l# 4 

5< 59 -260lg 163+H 145111 100±93° gglg 
Metal-Rich GCs ([Fe/H]> -0.905) 
~~0< 121 -25011^ 1541a 1 12ll? 191±H 

0-1 50 -2581J 4 1761?? 107l| 3151^ 19511 
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a Mean velocity. 

b Velocity dispersion. 

c Rotation-corrected velocity dispersion. 

d Rotation velocity derived from the rotation curve in the v p — vmzi 
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vs. X diagram. 

c Rotation velocity derived from the velocity histogram. 
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Table 2. Kinematics of the M31 Planetary Nebulae System 
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Fig. 1. — Comparison of the radial velocities for M31 GCs measured in this study with those 
in the literature (RBC2). The solid line indicates the linear least-squares fit, and the dashed 
line denotes the one-to-one relation. We labeled the ID of GCs in RBC2 for the GCs showing 
a velocity difference of over 110 km s _1 between this study and other studies. 
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Fig. 2. — Errors of the measured velocities of the M31 GCs versus T\ magnitude (a) and 
versus R, a galactocentric radial distance corrected for the inclination of M31 (b). In panel 
(a), filled circles for 136 metal-poor GCs derived in this study, open circles for 63 metal-rich 
GCs derived in this study, plusses for 174 metal-poor GCs derived in the previous studies, 
and crosses for 58 metal-rich GCs derived in the previous studies. In panel (b), filled circles 
for 46 bright GCs derived in this study, open circles for 99 faint GCs derived in this study, 
plusses for 155 bright GCs derived in the previous studies, and crosses for 104 faint GCs 
derived in the previous studies. 



-39- 



40 



e 

C\2 






3 
Z 



30 - 



20 - 



10 - 






1000-800 -600 -400 -200 

v [km s" 1 ] 







200 



Fig. 3. — Radial velocity histogram for all 504 GCs in the master catalog (open histogram) 
compared with 211 GCs measured in this study (hatched histogram). The vertical dotted 
line indicates the mean value of the radial velocities, v^ = —285 km s" 1 , for all GCs. 



-40- 



43° 00' 



30 



42° 00' 



o 
o 
o 
cv 



o 



30 



o 



g 41 u 00' 



o 

CD 
Q 



30 



40° 00' 



39° 30' 



L 1 


1 1 ' ' 


1 1 


1 1 1 1 1 1 


1 1 


1 1 . 


• 1 _ 


— 


> 










— 


- 




X 








- 


- ■ 


X 




X 






- 


- 


X * 










- 


t 




X 
X 


X 






- 








X 












O 


* O 




X 


- 


— 


X 


S x„ 


Xx* 


X 




— 


X 


O 

X 

x x 
x x 

°o 


CD 

xo 

I 
O *. 

x o 

o 


00„8U x x o 





X 

X 


- 




o 


X 


x*$$ x 4. x $£*£" x >< 








- 


• 


X 


Hf^« x 






- 


- 


» • 




< ■; > x x|^x§8 






- 


. 




X 


x x x ^ ~*x £V x x 




X 






- ■ 




° °« 
o P 




o 

O 




- 








«xx <^x 









- 






x ° X *■ 






- 






X 




• . x 






„ 




O , 


Ox x x ■ 








- 






* x x 
x x 






- 


- 






X 






- 


— 






X 






— 


- 






X 






- 


- 






X 






- 


- 










x 


x~ 


_ 






♦ 






_ 


r i 


, ' 1 , , 


i i 


1 1 1 II 1 1 1 


1 1 


1 . 


1 



00 h 50 m 



45 m 40 m 

Right Ascension (J2000) 



35 



m 



Fig. 4. — Positions of 211 GCs measured in this study (open circles) and of 293 GCs measured 
in the literature (crosses) excluding 61 GCs overlapping with our survey, overlaid on a 4° x 4° 
optical image of M31 from the Digitized Sky Survey. 
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the solid line indicates the sum of the two Gaussian fits, (b) Luminosity function for 483 
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magnitude limit to make subsamples of the bright and faint GCs, and the vertical dotted 
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Fig. 7. — Spatial distribution of all GCs with measured velocity (a), and radial velocities 
for M31 GCs versus the projected radii along the major axis (b) and along the minor axis 
(c). The receding GCs are plotted by open symbols, while the approaching GCs by filled 
symbols. The symbol size is relatively proportional to the velocity deviation. The large 
ellipse represents the optical extent of M31 based on the standard diam eter measured at a 
level of 25 mag arcsec -2 and ellipticity from iKarachentsev et al. I (120041 ) . Solid lines in (b) 
and (c) indicate the linear least-squares fits. 
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Fig. 8. — Radial velocities for M31 GCs as a function of projected distance along the major 
axis (left panels) and velocity histogram (right panels). Large open squares indicate the 
mean radial velocity of the GCs in a distance bin of 10 ' along the major axis. The vertical 
errorbar denotes the velocity dispersion of the GCs in the distance bin. In the velocity 
histogram, the sum of individual Gaussian fits is shown by a solid line, and the individual 
fits by dashed lines. The rotation fits for the GCs over the whole region are in (a) and (b), 
those in the range \Y\ < 1 kpc in (c) and (d), those in the range 1 < \Y\ < 3 kpc in (e) and 
(f), those in the range 3 < \Y\ < 5 kpc in (g) and (h), and those in the range 5 kpc < \Y\ in 
(i) and (j). Vertical dashed lines indicate the boundaries for the fit of a solid-body rotation. 
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Fig. 9. — Spatial distribution of the friendless GCs (a), and radial velocities versus the 
projected radii along the major axis (b), and along the minor axis (c). The receding GCs are 
plotted by filled circles, while the approaching GCs by crosses. The approaching friendless 
GCs are indicated by open squares, while the receding friendless GCs by open circles. The 
largest ellipse represents the optical extent of M31 based on the standard diameter measured 
at a level of 25 mag arcsec -2 and ellipticity. Satellite galaxies of M31 are indicated by small 
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Fig. 10. — Same as Fig. [91 but for the friendless GCs only. The symbol size in panel (a) is 
relatively proportional to the velocity deviation. 
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Fig. 11. — T\ magnitudes versus (C — T\) colors (a) and versus metallicities ([Fe/H]) (b), 
with histograms for T\ magnitudes (c), (C — T\) colors (d), and metallicities (e) for the 
friendless GCs. Normal GCs are indicated by small dots. The friendless GCs with smaller 
velocities than the rotation curve shown in Fig. [7J (b), are represented by open rectangles, 
and the friendless GCs with larger velocities than the rotation curve shown in Fig. [7] (b) by 
filled circles. The histogram for the friendless GCs is cross-hatched. 
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Fig. 12. — Spatial distribution of the metal-poor (filled circles) and metal-rich (crosses) GCs 
(a), and the histogram along the major axis (b), and along the minor axis (c). The largest 
ellipse represents the optical extent of M31 based on the standard diameter measured at 
a level of 25 mag arcsec -2 and ellipticity. Satellite galaxies of M31 are indicated by small 
ellipses. The histograms for the metal-poor and metal-rich GCs are indicated by open and 
hatched histograms, respectively. 






>- 



o 



>- 



-40 - 



100 



-49- 



— i 1 1 r~ 



- 1 1 1 1 - 



- 1 1 1 1— 



— 1 1 1 l~~ 




_l I I l_ 



' ' I L 



J I I l_ 



_l I I l_ 



50 







-50 



X [arcmin] 



-100 



Fig. 13. — (a) Number density contours for all the metallicity subsamples (solid lines for the 
metal-poor GCs and dashed lines for the metal-rich GCs), and (b) those for the the metallic- 
ity subsamples with T\ < 17. The contour levels are 0.01, 0.02, 0.03, 0.04 (GCs/arcmin 2 ) in 
(a) and 0.02, 0.04, 0.06, 0.09 (GCs/arcmin 2 ) in (b). The largest ellipse represents the optical 
extent of M31 based on the standard diameter measured at a level of 25 mag arcsec -2 and 
ellipticity. Satellite galaxies of M31 are indicated by small ellipses. 
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Fig. 15. — Radial velocities (dots) versus the projected distances along the major axis (X) 
and the velocity histogram for the MP (a,b) and MR (c,d) GCs. Large open squares indicate 
the mean radial velocity of the GCs in a distance bin of 10 ' along the major axis. The vertical 
errorbar denotes the velocity dispersion of the GCs in the distance bin. The contours indicate 
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Fig. 16. — Same as Fig. [HI but for the metal-poor GCs. 
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Fig. 17. — Same as Fig. [BJ but for the metal- rich GCs. 
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Fig. 18. — Same as Fig. [121 but for the bright (filled circles and open histogram) and faint 
(crosses and hatched histogram) GCs. 
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Fig. 20. — Same as Fig. [15], but for the bright (a,b) and faint (c,d) GCs 
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Fig. 21. — Spatial distribution of X-ray emitting GCs (a), and radial velocities versus the 
projected radii along the major axis (b) and along the minor axis (c). We plot 39 genuine 
GCs with X-ray detection and measured velocity using variable symbol sizes according to 
their velocity deviations from the systemic velocity of M31 (receding GCs with open circles 
and approaching GCs with filled rectangles). One GC and 16 GC candidates with X-ray 
detection but no measured velocity are plotted by crosses, and all 504 GCs by dots. The 
largest ellipse represents the optical extent of M31 based on the standard diameter measured 
at a level of 25 mag arcsec -2 and ellipticity. Satellite galaxies of M31 are indicated by thick 
pluses with small ellipses. Receding and approaching X-ray emitting GCs with measured 
velocity are plotted by open circles and rectangles, respectively, in (b) and (c). 
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Fig. 23. — Rotation curves (a) and velocity dispersion profiles (b) derived using the data 
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Fig. 25. — Mean metallicity versus the projected galactocentric radi us for the metal-poor 



GCs (starlets) in comparison with the red giant branch stars given by iKalirai et al. I ((2006) 
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(circles). The vertical error bars represent the mean errors (solid line) and dispersion (dashed 
line) of the mean metallicity and the horizontal error bars represent the radial coverage of 
the data. 
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Fig. 2 6. — Same as Figure 8, but for the PNe based on the data given by iMerrett et al. 
<J2006h . 



